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Abstract. We give a strong evidence that noncrystalline materials such as quasicrystals 
or incommensurate solids are not exceptions but rather are generic in some regions of a 
phase space. We show this by constructing classical lattice gas models with translation- 
invariant, finite-range interactions and with a unique quasiperiodic ground state which is 
stable against small perturbations of two-body potentials. More generally, we provide a 
criterion for stability of nonperiodic ground states. 
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1 Introduction 



One of the important problems in physics, the so-called crystal problem, is to understand 
why matter is crystalline at low temperatures [|T], |^, |3], |, It is traditionally 

assumed (but has never been proved) that at zero temperature, minimization of the free 
energy of a system of many interacting particles can only be obtained by their periodic 
arrangements (a perfect crystal), which at nonzero temperature are disrupted by defects 
due to entropy. Recently, however, there has been a growing evidence that this basic 
phenomenon, the crystalline symmetry of low temperature matter, has exceptions; in 
particular incommensurate solids [§] and, more recently, quasicrystals ||. It is very 
important to find out how generic these examples are. In other words, is nonperiodic 
order present in these systems stable against small perturbations of interactions between 
particles? 

The problem of stability of quasiperiodic structures was studied recently in continuum 
models of particles interacting through a well, Lennard- Jones, and other potentials |]T0 



TT| , [T2| , |13| , |14| , [L5| . However, no final conclusion was reached. After all, one has to compare 
chosen quasiperiodic structure with all possible arrangements of particles in the space, a 
really formidable task. 

Here we will present two classical lattice gas (toy) models with a unique stable non- 
periodic ground state. More precisely, every site of a square lattice can be occupied by 
one of several different particles. The particles interact through two-body finite-range 
translation-invariant potentials (we chose chemical potentials of all particles to be zero). 
Our models have only nonperiodic ground-state configurations (infinite-lattice configura- 
tions minimizing the energy density of the system) . We will prove that if one perturbs our 
models by sufficiently small chemical potentials or two-body translation-invariant inter- 
actions, their ground-state configurations do not change. It means that there is an open 
set in a space of two-body interactions without periodic ground-state configurations. This 
constitutes a first generic counterexample to the crystal problem. 

In Section 2, we discuss general classical lattice gas models with a unique nonperiodic 
ground state. We introduce a criterion, the so-called strict boundary condition, for nonpe- 
riodic ground- state configurations to be stable against small perturbations of finite-range 
potentials. In Section 3, we describe main features of Robinson's nonperiodic tilings of 
the plane, construct a classical lattice gas model based on them, and show why its unique 
nonperiodic ground state is not stable against small perturbations of chemical potentials. 
In Section 4, we present modifications of Robinson's tilings which allows us to construct 
models with a unique nonperiodic ground state satisfying our stability criteria. Section 5 
contains a short discussion. 



2 Classical Lattice Gas Models and Nonperiodic - 
Ground States 

A classical lattice gas model is a system in which every site of a lattice Z d can be occupied 
by one of n different particles. An infinite-lattice configuration is an assignment of parti- 
cles to lattice sites, i.e., an element of Q = {1, ...,n} z . If X e Q and A C Z d , then we 
denote by X (A) a restriction of X to A. Particles at lattice sites a and b interact through 
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a two-body translation- invariant potential $(a — b), which is a function on {1, n}^ a ' h ^ 
- the space of all possible assignments of particles to lattice sites a and b, and we assume 
that $(a — b) = if dist(a, b) > r. The Hamiltonian in a bounded region A can be then 
written as H\ = J2a,beA ^( a ~ b). 

Y is a local excitation of X, Y ~ X, Y, X 6 , if there exists a bounded A C Z d 
such that X = Y outside A. 

For Y ~ X the relative Hamiltonian is defined by 

H(Y,X)= £ ($(a-b)(F)-$(a-b)(X)). 

{a,b}nA^0 

X 6 is a ground-state configuration of H if 

#(Y, X) > /or any Y ~ X. 

That is, we cannot lower the energy of a ground-state configuration by changing it 
locally. 

The energy density e(X) of a configuration X is 

e(X) = hminf^P, 
v ' A->Z<* |A| 

where |A| is the number of lattice sites in A. It can be shown that any ground-state 
configuration has the minimal energy density (for a proof see Appendix A). It means 
that local conditions present in the definition of a ground-state configuration force global 
minimization of the energy density. 

Equilibrium behavior of a system of many interacting particles can be described by 
a grand-canonical ensemble which tells us what are probabilities of finding particles at 
given positions. An infinite- volume limit of this ensemble is called a translation-invariant 
Gibbs state or equilibrium state. Mathematically speaking, it is a measure on the space, 
Q, of all configurations (for precise definitions see Appendix B). Then by a ground state 
we mean a limit of a nonzero- temperature equilibrium state as temperature approaches 
zero, with other variables, such as chemical potentials, kept fixed. A ground state is, 
therefore, a translation-invariant probability measure on Q. It can be shown that the set 
of all ground- state configurations have probability one with respect to a ground-state 
measure. 

If a system has a unique periodic ground-state configuration and its translations, then 
a unique ground-state measure assigns an equal probability to all these translations. For 
example, the Ising antiferromagnet has two alternating ground-state configurations but 
only one ground-state measure which assigns probability 1/2 to both of them. 

Generally, a probability ground-state measure \i gives equal weights to all ground-state 
configurations and can be obtained as a limit of averaging over a given ground-state config- 
uration X and its translations r a X by lattice vectors a G Z d : /i = lim A ^ Z d 4r X) a eA ^( r aX), 
where 5(r a X) is a probability measure assigning probability 1 to r a X. 

For any potential, the set of ground-state configurations is nonempty but it may not 
contain any periodic configurations |1| [L7], [18], [19], |2(| [21], ^3] . We restrict ourselves 



to systems in which, although all ground-state configurations are nonperiodic, there is 
a unique translation-invariant probability measure supported by them. Therefore, such 
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measure is inevitably a zero-temperature limit of translation-invariant Gibbs states (equi- 
librium states) and hence it is a ground state of a given model. Uniqueness of a ground- 
state measure is equivalent to the statement that a uniformly defined frequency of any 
finite arrangement of particles is the same in all ground-state configurations (such mea- 
sures are called uniquely ergodic measures and their supports are called uniquely ergodic 
sets) As a support of a ground-state measure we will always choose a minimal set 
(if any arrangement appears in a ground-state configuration, it appears with a positive 
density), i.e., we exclude ground-state configurations with broken bonds (such as interface 
ground-state configurations of the ferromagnetic Ising model) . More precisely, to find the 
frequency of a finite arrangement in a given configuration we first count the number of 
times it appears in a box of size I which is centered at the origin of the lattice, divide 
it by l d , and then take the limit I — > oo. If the convergence is uniform with respect to 
the position of the boxes, then we say that the configuration has a uniformly defined 
frequency of this arrangement. 

Let us emphasize again that in our models the frequency of any finite particle ar- 
rangement is the same in all ground-state configurations. Yet, it is not true that all 
ground-state configurations are lattice translations of a single nonperiodic ground-state 
configuration. That is why, in order to deal with all ground-state configurations simul- 
taneously, we have to introduce a ground-state measure. Such a situation is present in 
all tiling models of quasicrystals; for example in the Penrose tilings, where instead of 
ground-state configurations we have perfect tilings. 

To get all ground-state configurations we have to close (in a topology described in 
Appendix B) a set of all translations of any ground-state configuration (the closure does 
not depend on the choice of a configuration). In this way we obtain uncountably many 
ground-state configurations. In fact, probability of any single configuration should be zero; 
otherwise, the measure of a countable union of its translations would be infinite. If there 
were countably many ground-state configurations (there are countably many translations), 
then the measure of the set of all ground-state configurations would be zero, and this is 
impossible. 

We will introduce now a condition which is equivalent to the stability of nonperiodic 
ground states. It generalizes the so- called Peierls condition for models without 
periodic ground-state configurations. 

For clarity of presentation we assume that our models are two- dimensional. We 
also assume that our models have a unique ground state, supported by ground-state 
configurations for which all interactions attain simultaneously their minima (we choose 
them to be equal to zero). As it was explained before, we do not have to consider 
ground-state configurations with broken bonds. Therefore, if Y is not a ground-state 
configuration, it contains pairs of particles with nonminimal energies (we choose them 
to be equal to 1), the so-called broken bonds. Denote by B{Y) the number of broken 
bonds in Y . Thus H(Y, X) = B(Y) if Y ~ X and X is a ground-state configuration. 

Condition The strict boundary condition for local excitations. 

Let X be a ground-state configuration and Y a local excitation of X; Y ~ X. Let 
n ar (Y,X) denote the difference of the number of appearances of an arrangement ar (a 
particle or a pair of particles) in Y and the number of its appearances in X. 

We say that a model satisfies the strict boundary condition for local excitations and 
an arrangement ar if there exists a C ar > such that for every ground-state configuration 
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X and every Y ~ X 

\n ar (Y,X)\<C ar B(Y). (1) 

Theorem 1 A unique ground state of a finite-range Hamiltonian, H , is stable against 
small perturbations of chemical potentials and two-body interactions of range smaller than 
r if and only if the strict boundary condition is satisfied for particles and pairs of particles 
at a distant smaller than r . 

Proof: Assume first that the strict boundary condition does not hold for a ground-state 
configuration X and an arrangement ar (a particle or a pair of particles at a distance 
smaller than r). Hence, for any C > there exists Yq ~ X such that n ar (Yc, X) > 
CB{Yc) or n ar (Yc, X) < —CB(Yc)- We assume that the first case holds; the second one 
can be treated in an analogous way. We introduce a perturbation, which assigns a 
negative energy, E = — 1/C, to the arrangement ar. This perturbation is small if C is 
large. Let H" = H + H' , where H' = Ea *a- Then 

H"(Y C , X) = B{Y C ) - n ar (Y c , X)/C < 0. (2) 

Thus, X is not a ground-state configuration for the perturbed Hamiltonian H" . 

Assume now that \n ar (Y, X)\ < C ar B(Y) for every arrangement ar (a particle or a pair 
of particles at a distance smaller than r) and every local excitation Y ~ X of a ground- 
state configuration X. Let \P be a potential of a range smaller than r and H" = H + H' 
be a perturbed Hamiltonian, where H' = Ea^a- Then 

H"(Y,X) = B(Y)+J2nar(Y,X)V(ar) (3) 

ar 

>B(Y)-B(Y)-£C ar \V(ar)\>0 (4) 

ar 

if 

£C or |*(ar)|<l. (5) 

ar 

Hence, X is a ground-state configuration for if" for every sufficiently small perturbation. 

So far we have proved that ground-state configurations of H (without broken bonds) 
are ground-state configurations of H" . Now we will show that the set of ground-state 
configurations of H" is a uniquely ergodic set. It will then follow that the ground state 
of H" must be equal to that of H . Let fj," be an ergodic measure different from the 
ground-state measure of H. It has therefore a nonzero density, u>, of broken bonds. Let 
Z be in the support of //' and have the density uo of broken bonds. We will now show 
that Z is not a ground-state configuration of H" . 

We construct two configurations: W ~ Z such that W = X on a square A of size / 
and W = Z outside A, and Y ~ X such that Y = Z on A and Y = X outside A. 

Since Y is an excitation of X, then 

K r (F,X)| < C ar B(Y) < C ar (ul 2 + o(l 2 )). (6) 
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When we change Z to X on A we decrease the number of broken bonds at least by 
ul 2 — o{l 2 ) and change the number of appearances of ar at most by C ar {ujl 2 + o{l 2 )) which 
follows from (7). Hence, 

H"(W, Z) < -(ul 2 - oil 2 )) + £ \*(ar)\(C ar (ul 2 + o(/ 2 )). (7) 

ar 

Thus, 

H"(W, Z) < (8) 

if (6) holds and Z is big enough. It follows that Z is not a ground-state configuration of 
H". □ 



3 Robinson's Tilings and Unstable Ground States 

To construct our classical lattice gas models, we will use Robinson's nonperiodic tilings 
of the plane [p8 , [29|j . He designed 56 square-like tiles such that using an infinite number 
of their copies one can tile the plane only in a nonperiodic fashion. We describe now the 



Robinson tiles; we follow p8| closely. There are five basic tiles represented symbolically 
in Fig.l; all other tiles can be obtained from them by rotations and reflections. The first 
tile on the left is called a cross; the other tiles are called arms. Every tile is also furnished 
with one of the four parity markings shown in Fig. 2. Crosses can be combined with the 
parity marking at the lower left in Fig.2. Vertical arms (the direction of long arrows) can 
be combined with the marking at the lower right and horizontal arms with the marking at 
the upper left. All tiles can be combined with the remaining marking. Two tiles "match" 
if arrow heads meet arrow tails separately for the parity markings and the markings of 
crosses and arms. 

We will now describe main features of Robinson's nonperiodic tilings. Observe first 
that, in any tiling, the centers of the tiles form the square lattice Z 2 . Therefore, any tiling 
can be described by an assignment of tiles to lattice sites, a lattice configuration. It is also 
easy to see that if the plane is tiled with tiles with the above parity markings, then these 
must alternate both horizontally and vertically in the manner shown in Fig.2. Hence, 
every odd-odd position on the Z 2 lattice (if columns and rows are suitably numbered) is 
occupied by crosses with the lower left parity marking. They form a periodic configuration 
with period 4 as shown in Fig.3 (crosses are denoted there by |_, where line segments denote 
double arrows). Then in the center of each "square" matching rules force a cross such 
that the previous pattern reproduces but this time with period 8. 

Continuing this procedure (the so-called principle of expanding squares) infinitely 
many times, we obtain a nonperiodic configuration. It has built in periodic configurations 
of period 2 n , n > 2, with tiles on sublattices of Z 2 as is shown in Fig. 4. Each side of 
every square consists of two lines of arms with long double arrows at the left of Fig.l 
(called double arms) meeting half-way through double short arrows of arms at the top of 
Fig.l. It is exactly the arms with asymmetrically located double short arrows which force 
crosses to be present at the center of each square. Crosses of different squares at the same 
sublattice are joined by lines of arms with single long arrows at the right of Fig.l (called 
single arms). For a complete construction of Robinson's tilings and proofs, we refer the 
reader to the original article of Robinson [28|. 
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Now, using Robinson's tiles we will construct a lattice gas model WB, O, 113, IT9L pOL ETl 



22| . |23fl . Every site on a square lattice can be occupied by one of the particles corresponding 
to the tiles. Two nearest-neighbor particles which do not match as tiles have a positive 
interaction energy, say 1; otherwise, the energy is zero. We obtain a lattice gas model 
with nearest-neighbor translation-invariant interactions. 

Any periodic configuration with period p corresponds to a periodic configuration of 
tiles. Therefore, it has, in any square of p 2 lattice sites, at least one pair of nearest 
neighbors with energy equal to 1. Hence, any periodic configuration has a nonzero energy 
density. On the other hand, nonperiodic configurations corresponding to tilings have 
energy density equal to zero. It follows that the model does not have periodic ground- 
state configurations. 

It has been proven [jnj that the unique ground state of this model is not stable against 
small perturbations of chemical potentials. To see this let us introduce a small negative 
chemical potential h for single arms. Notice that every ground- state configuration con- 
tains arbitrarily long sequences of double arms. When we change such a sequence of length 
I into a sequence of single arms, we lower the energy by Ih along the sequence, and only 
at two endpoints of the sequence the energy increases by 1. Hence, no matter how small 
h is, one can always change locally a ground-state configuration of an unperturbed model 
to lower its energy. It is easy to see that the strict boundary condition is not satisfied. 

In fact, we showed that frequencies of single and double arms change continuously with 
respect to their chemical potentials. This behavior of continous change in stoichiometries 
at zero temperature is not present in previously known models and real alloys. 



4 Stable Quasiperiodic Ground States 

We will present in a moment models with a unique nonperiodic ground state satisfying 
our stability criteria. This will be a modification of the model based on Robinson's tiles 
such that a sufficiently small chemical potential for one kind of arm and a small two-body 
perturbation do not change its unique ground state. 

We cannot simply get rid of the double arms because it is exactly the double arms 
which create squares whose central crosses are forced. The absence of double arms will 
destroy the principle of expanding squares and we will no longer have a unique nonperiodic 
ground state. 

However, let us recall that the forcing is done by four arms located exactly at the 
middle point between crosses forming a square; these are the arms at the top of Fig.l. 
The main idea of our construction is to force similar arms in the middle of every side 
of every square, as in the original Robinson tiling, by using periodic sequences of new 
markings instead of translation- invariant sequences of single and double arms. We will 
use sequences of alternating A, B, and C markings propagating from crosses. Whenever 
two B or C markings meet at the middle of a side of a square, then our special arms are 
forced there. 

Observe that we do not have long sequences of different arms anymore, just three 
translations of lines ...ABC ABC... It is impossible to distinguish between these transla- 
tions by chemical potentials and also by two-body interactions because of the symmetry of 
the tiles and the fact that 2 n "periods" of Robinson's tilings and period 3 of the sequences 
are relatively prime. 
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MODEL I 



We will describe first our tiles. These are modified crosses and arms with the parity- 
markings as before. We have also two additional levels of letter markings. Each side of 
each tile is marked by A, B, or C and by T, O, M, E, or K (not all choices are present 
as explained below). All possible markings of crosses are shown in Fig. 5. Next we allow 
all rotations of crosses but we do not allow any reflections. This feature makes our model 
somewhat different from the original Robinson one. We have 15 types of arrows present 
in arms: all combinations of A — > B ^ C, C ^ A and T — > O, O — > M, M — > E, 
E — > K, K — > T. Line segments present in arms are marked A — A, B — B, or C — C 
while B — B and C — C segments are also marked T — O or O — T as is shown in Fig. 6. 
Again, we allow rotations but not reflections. 

Two tiles match, if in addition to previous requirements, letter markings match sepa- 
rately on both levels. We also assume that line segments A — A can be put next to any 
letter marking of the second level (this is not a tiling condition). 

Our matching rules can be translated as before into a lattice gas model with nearest- 
neighbor interactions. In addition, we have to take care of the fact that letter markings 
of the second level which are opposite nearest neighbors of A — A line segments are not 
forced to be correlated. We introduce a two-body interaction of range 2 between two 
arrows facing one another such that the energy is zero if one of them has T marking and 
the other one is marked by O; otherwise, the energy is 1. 

We will show that the unique ground-state measure of the above model is of the Robin- 
son type. More precisely, looking at crosses and ignoring A arrows on every sublattice 
2 n Z 2 with n odd and C arrows for n even, we see the self-similar structure of expanding 
squares shown in Fig. 4. 

Proposition 1 The unique ground-state measure of Model I is of the Robinson type as 
described above. 

Proof: First, let us recall that the parity markings force crosses to occupy every site of a 
2Z 2 sublattice. Letter markings of the second level allow only crosses which are rotations 
of the cross at the upper left in Fig. 5. 

Arms with the upper left and lower right parity marking force relative orientations of 
these crosses to be such that squares of length 2 are created as is shown in Fig. 3 (arrows 
of crosses with A markings and arms are not drawn and remaining arrows are represented 
by line segments). The above arms with C — C line segments force crosses to appear in 
the center of every square just created. New crosses (four rotations of the cross at the 
upper right in Fig. 5) form a 4Z 2 sublattice. 

Once we have crosses on a sublattice, lines of arrows propagate from every cross and 
meet half-way through line segments of arms. For the crosses on a 2 n Z 2 sublattice with n 
odd, these are alternating C — C and A — A line segments and alternating B — B and A — A 
segments for n even. These segments force crosses to have relative orientations on every 
sublattice as is shown in Fig.3. C — C(B — B) segments force crosses in centers of squares 
just created. Crosses on successive sublattices have to be chosen in clockwise order from 
Fig.5. The upper left corner of every square has to be occupied by an appropriate cross 
from Fig.5 and other corners by clockwise rotations. An infinite nonperiodic configuration 
is inductively created. □ 
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Proposition 2 Model I satisfies the strict boundary condition for particles and pairs of 
particles, if we do not distinguish between their different rotations and reflections. 

Proof: Let a broken bond be a unit segment on the dual lattice separating two nearest- 
neighbor particles with a positive interaction energy (a common side of two nearest- 
neighbor tiles which do not match) or separating two arrows at a distance 2 with the 
positive energy. Let us then divide an excitation Y of a ground-state configuration X 
into connected components without broken bonds (two lattice sites are connected if they 
are nearest neighbors) such that on every component crosses having relative orientations 
shown in Fig. 3 form a 2Z 2 sublattice. This is achieved by paths on the lattice dual to Z 2 
with lengths not bigger than 16 and joining broken bonds. 

When one considers only 2Z 2 sublattices and lines of arms joining their sites, the 
absolute value of the difference between the number of appearances of a given particle in 
Y and the number of its appearances in X is not bigger than 16 times the number of the 
above paths. 

Now we decompose every connected component into components with crosses that 
form 4Z 2 sublattices. This is achieved by paths on the dual lattice with lengths not 
bigger than 32 and again joining broken bonds. When one considers 4Z 2 sublattices and 
lines of arms joining their sites, the additional difference of the number of appearances 
of our particle is not bigger than 16 times the number of new paths plus 16/2 times the 
number of old paths. 

At the n — th step we divide every connected component constructed so far into 
components with crosses forming 2 n Z 2 sublattices and we use paths with lengths not 
bigger than 2 n+3 . Again, the additional difference of the number of appearances of our 
particle is bounded by 16 times the number of new paths plus 1/2 of the difference at the 
(n — l)th step. 

We repeat this process m times, where m is the smallest number such that diam{&, Y(a)- 
7^ A (a)} < 2 m . Now, broken bonds and paths can be regarded as vertices and edges of 
a planar graph. It follows from Euler's formula that the total number of all paths is 
bounded by 3 times the number of broken bonds. This shows that \n(Y,X)\ < 9QB(Y) 
for a given particle. 

We will now discuss two-particle arrangements. For a pair of particles at a distance 
D, at least one particle is located on a sublattice 2 fe Z 2 with 2 k < 2D or on a line joining 
sites of this sublattice. To get a bound for a pair of particles we must take into account 
the effect of shifting one line of arms joining two sites of a sublattice with respect to the 
other parallel line. This can be done by multiplying the previous bound by 2D. □ 

Let us now perturb our model. We will consider only covariant interactions. That is, 
a chemical potential of a particle-tile or a two-body interaction between a pair of particles 
should be the same for any of their rotations and reflections. Combining Theorem 1 and 
Proposition 2 we get 

Theorem 2 The unique nonperiodic ground state of Model I is the unique ground state 
for any sufficiently small chemical potential and two-body covariant perturbation. 

MODEL II 

Let us now describe tiles of our second model. We have modified crosses with A and C 
arrows shown in Fig. 7, arms with A — > B, B — > C, and C — > A arrows and A — A, B — B, 
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and C — C line segments shown in Fig. 8 and 9. Arms in Fig.8 are combined with the 
upper right parity marking. Arms in Fig.9 can be combined with the lower right or upper 
left parity marking. Crosses are combined with the parity markings as before. We allow 
all rotations and reflections of our new tiles as in the original Robinson model. 

We do not have anymore the previous second level of letter markings. However, B — B 
and C — C line segments are not located in the middle of arms anymore so they force 
crosses to create squares. Arms with central A — A segments from Fig.8 cannot prevent 
these squares from going out of phase; a fault shown in Fig. 10 can be created. Such faults 
are also present in some of the original Robinson tilings. One may prove that any tiling 
has at most one fault (follow the proof of the analogous theorem in [28], look also at the 



proof of Proposition 3). Therefore, as in the previous model, the unique ground state of 
Model II is of the Robinson type. 

Proposition 3 Model II satisfies the strict boundary condition for particles and pairs of 
nearest-neighbor particles, if we do not distinguish between their different rotations and 
reflections. 

Proof: Observe first that arms with the lower right and the upper left parity marking 
have asymmetric A — A line segments. Therefore, there are no faults on 2Z 2 sublattices. 
The first step is the same as in the proof of Proposition 2. 

In the second step we decompose every connected component without broken bonds 
into components with crosses that form 4Z 2 sublattices but not necessarily having relative 
orientations as in Fig. 3; faults may be present. This is achieved by paths on the dual lattice 
with lengths not bigger than 32 and joining broken bonds. 

We may have an arbitrarily long faults perpendicular to A — A line segments without 
any broken bonds except at the endpoints of the fault. The presence of two arms on the 
fault with line segments parallel to the fault forces a cross (or a broken bond) between 
them and then a broken bond is forced at a distance not bigger than 16 from the cross. 
Arms with B — B or C — C line segments perpendicular to the fault also force a cross on 
the fault or a broken bond at a distance not bigger than 16. Now, when one considers 
only 4Z 2 sublattices and lines of arms joining their sites, then the absolute value of the 
difference between the number of appearances of a given particle in an excitation Y and 
the number of its appearances in a ground-state configuration X is not bigger than 16 
times the number of paths in the present decomposition plus 16/2 times the number of 
paths in the first step (there is no change of the number of particles and pairs of nearest 
neighbors along a fault with A — A line segments perpendicular to it). 

Now we decompose every connected component into components with crosses that 
form 8Z 2 sublattices. These crosses go out of phase along previous faults and they may 
also create they own faults. We repeat this process m times, as in the proof of Proposition 
2, to get \n(Y, X)\ < 965(F) for particles and pairs of nearest-neighbor particles. □ 

Theorem 3 The unique nonperiodic ground state of Model II is the unique ground state 
for any sufficiently small chemical potential and nearest-neighbor covariant perturbation. 

Let us notice that particles of our models have shapes and interactions between them 
are covariant but not symmetric: an interaction between nearest neighbors RP is the 
same as between but not necessarily the same as between PR. There is, however, 



10 



a construction due to Radin |L7j who enlarges a family of possible particles and makes 
them shapeless, and interactions between them fully symmetric (dependent only on the 
distance between given particles). He constructs in this way a classical lattice gas model 
with a finite-range fully symmetric interactions and with a unique nonperiodic ground 
state. His method requires the original family of particles-tiles to be both rotation and 
reflection-invariant. Therefore, it can be applied only to our second model. We use it to 
construct a model with shapeless particles and with a unique nonperiodic ground state 
which is stable in a space of fully symmetric interactions. 

Theorem 4 The unique ground state of the symmetric version of Model II is the unique 
ground state for any sufficiently small chemical potential and symmetric perturbation of 
range smaller than three lattice spacings. 



5 Conclusions 

We have constructed two finite-range classical lattice gas models with a unique quasiperi- 
odic ground state which is stable against small perturbations of interactions. 

The family of particles of our first model is rotation-invariant but not reflection- 
invariant. Our quasiperiodic ground state is stable in the space of covariant two-body 
interactions. 

In the second model, particles are shapeless and interaction between them depends 
only on the distance between them. The unique quasiperiodic ground state is stable in 
the space of fully symmetric interactions of range smaller than three lattice spacings. 

In the construction of our models we used the presence of order in nonperiodic ground- 
state configurations. Although nonperiodic, they are not, however, completely chaotic. 
They exhibit long-range positional order in the sense that arrangements of particles at 
distant regions are correlated. In fact, our ground- state configurations possess highly 
ordered structures: if a certain fraction of particles is ignored, the rest of a ground state 
configuration is periodic; the smaller the fraction, the larger the period. It would be 
interesting to see if other, more disordered, ground states are stable. Modifying the 
principle of expanding squares responsible for nonperiodicity of Robinson's tilings, many 
models without periodic ground states were recently constructed [[30], |3|, |3l| . 



Let us point out that the Robinson tilings share with the Penrose tilings many essential 
features like forced nonperiodicity, self- similarity, uniqueness of a structure for given 
matching rules. In fact, using the so-called Ammann bars, the Penrose tilings can be 



also represented by configurations on a square lattice [p9|| . We hope that investigating 
Robinson's tilings one may provide some information about behavior of quasicrystals 
associated with the Penrose tilings. Let us also note that models based on the Penrose 



tilings and with particles on a square lattice were recently investigated in Ref. ||32 |. 

Finally, let us say few words about positive temperatures, where to construct stable 
equilibrium phases of a system of many interacting particles one has to minimize not the 
energy but the free energy of the system taking into account disruptions due to entropy. 
For lattice gas models with a finite number of periodic ground-state configurations satis- 
fying the so-called Peierls condition, there exists a complete theory due to Pirogov and 
Y.Sinai |2^, ^B, 27J. It tells us that for every Hamiltonian with the above properties, a 
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low-temperature phase diagram is a small deformation of a zero-temperature phase di- 
agram. We would like to generalize the Pirogov-Sinai theory to be able to apply it to 
lattice models without any periodic ground-state configurations. One of the important 
problems is to construct a quasiperiodic equilibrium phase for a finite-range interaction 
- a microscopic model of a quasicrystal. We conjecture the existence of quasiperiodic 
equilibrium phases in a modified three-dimensional version of our second model. 

6 Appendix A 

Theorem 5 Let X be a ground-state configuration of a finite-range Hamiltonian H . Then 
X has the minimal energy density: e(X) < e(Y) for any FgO. 

Proof by contraposition: For clarity of presentation we consider two-dimensional mod- 
els. Assume that there exists a configuration Y with smaller energy density than that 
for a ground-state configuration X: e(Y) < e(X). Then there is a square, Sq(N), of N 2 
lattice sites and 5 > such that 

Hl q{N) (Y) < Hl q(N) (X) - N 2 5 (9) 

and 

8N^max We n Lambda \MW)\ < N 2 5. (10) 
oeA 

Now we construct an excitation Z ~ X as follows: Z(Sq(N)) = Y(Sq(N)) and Z = X 
outside Sq(N). It follows that H(Z,X) < so X is not a ground-state configuration. 

7 Appendix B 

We equip {1, n} with the discrete topology (every subset of {1, n} is both open and 
closed). Q is then compact in the product topology. Let ar be a finite arrangement of 
particles, ar e {1, ...,n} A for some bounded A, then 

C a A r = {X eVl: X(A) = ar} 

is a cylinder set based in A. Cylinder sets are open and they generate all open sets in Q 
(every open set is a union of some cylinders). It follows that a sequence of configurations 
X m converges to X, 

X m — > X as m — > oo, 

if for every bounded A e Z d there is m(A) such that for every m > m(A), X m (A) = X(A). 

Borel sets are generated from open sets by taking complements and countable unions. 

A probability measure, /i, on a configuration space, Q, is a function which assigns 
a number between and 1 to every Borel set such that = 1, /x(0) = 0, and 

oo 

1=1 

for Borel sets Ai and Aj if Ai n A,- = for all « 7^ j. 



12 



Acknowledgments. I have begun this research at Institut de Physique Theorique, Unite 
FYMA, Catholique Universite de Louvain-la-Neuve (it was supported there by Bourse de 
recherche UCL/FDS) and continued at Instituut voor Theoretische Fysica, Katholieke 
Universiteit Leuven. I would like to thank the Polish Committee for Scientific Research - 
KBN for a financial support under the grant 2P03A01511. 

References 

[1] P. W. Anderson, Basic Notions of Condensed Matter Physics (Benjamin, Menlo Park, 
1984), p.12. 

[2] S. G. Brush, Statistical Physics and the Atomic Theory of Matter, from Boyle and 
Newton to Landau and Onsager (Princeton University Press, Princeton, 1983), p. 277. 

[3] B. Simon, in Perspectives in Mathematics: Anniversary of Oberwolfach 1984 
(Birkhauser-Verlag, Bassel, 1984), p.442. 

[4] G. E. Uhlenbeck, in Fundamental Problems in Statistical Mechanics II, edited by E. 
G. D. Cohen (Wiley, New York, 1968), p. 16. 

[5] G. E. Uhlenbeck, in Statistical Mechanics; Foundations and Applications, edited by 
T. A. Bak (Benjamin, New York, 1967), p.581. 

[6] C. Radin, Low temperature and the origin of crystalline symmetry, Int. J. Mod. Phys. 
Bl: 1157 (1987). 

[7] C. Radin, Global order from local sources, Bull. Amer. Math. Soc. 25: 335 (1991). 

[8] S. Aubry, Devil's staircase and order without periodicity in classical condensed mat- 
ter, J. Physique 44: 147 (1983). 

[9] Aperiodicity and Order Vol 1 and 2, edited by M. V. Jaric (Academic Press, London, 
1987, 1989). 

[10] M. Widom, K. J. Strandburg, and R. H. Swendsen, Quasicrystal ground state, Phys. 
Rev. Lett. 58: 706 (1987). 

[11] S. Narasimhan and M. V. Jaric, Icosahedral quasiperiodic ground states? Phys. Rev. 
Lett. 62: 454 (1989). 

[12] A. P. Smith and D. A. Rabson, Comment on "Icosahedral quasiperiodic ground 
states?" Phys. Rev. Lett. 63: 2768 (1989). 

[13] S. Narasimhan and M. V. Jaric reply Phys. Rev. Lett. 63: 2769 (1989). 

[14] Z. Olami, Stable dense icosahedral quasicrystals, Phys. Rev. Lett. 65: 2559 (1990). 

[15] A. P. Smith, Stable one-component quasicrystals, Phys. Rev. B 43: 11635 (1991). 

[16] C. Radin, Tiling, periodicity, and crystals, J. Math. Phys. 26: 1342 (1985). 



13 



[17] C. Radin, Crystals and quasicrystals: a lattice gas model, Phys. Lett. 114A: 381 
(1986). 

[18] C. Radin, Crystals and quasicrystals: a continuum model, Commun. Math. Phys. 
105: 385 (1986). 

[19] J. Miekisz and C. Radin, The unstable chemical structure of the quasicrystalline 
alloys, Phys. Lett. 119A: 133 (1986). 

[20] J. Miekisz, Many phases in systems without periodic ground states, Commun. Math. 
Phys. 107: 577 (1986). 

[21] J. Miekisz, Classical lattice gas model with a unique nondegenerate but unstable pe- 
riodic ground state configuration, Commun. Math. Phys. Ill: 533 (1987). 

[22] J. Mi§kisz, A microscopic model with quasicrystalline properties, J. Stat. Phys. 58: 
1137 (1990). 

[23] C. Radin, Disordered ground states of classical lattice models, Rev. Math. Phys. 3: 
125 (1991). 

[24] J. C. Oxtoby, Ergodic sets, Bull. Amer. Math. Soc. 58: 116 (1952). 

[25] Pirogov and Y. G. Sinai, Phase diagrams of classical lattice systems I, Theor. Math. 
Phys. 25: 1185 (1975). 

[26] Pirogov and Y. G. Sinai, Phase diagrams of classical lattice systems II, Theor. Math. 
Phys. 26: 39 (1976). 

[27] Ya. G. Sinai, Theory of Phase Transitions: Rigorous Results (Pergamon Press, Ox- 
ford, 1982). 

[28] R. M. Robinson, Undecidability and nonperiodicity for tilings of the plane, Invent. 
Math. 12: 177 (1971). 

[29] B. Griinbaum and G. C. Shephard, Tilings and Patterns (Freeman, New York, 1986). 

[30] S. Mozes, Tilings, substitutions, and dynamical systems generated by them, J. 
dAnalyse Math. 53: 139 (1989). 

[31] B. Tsirelson, A robust nonperiodic tiling system, preprint of Tel Aviv University, 
School of Mathematical Sciences (1992). 

[32] W. Li, H. Park, and M. Widom, Phase diagram of a random tiling quasicrystal, J. 
Stat. Phys. 66: 1 (1992). 



14 



